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Длинные кости конечностей четвероногих как система рычагов, 
обеспечивающих статику и локомоцию тела, постоянно подвержены воз- 
действию различных по силе и характеру статических и динамических 
нагрузок., Реакция на эти нагрузки проявляется в перестройке архитек- 
туры, изменений формы и размеров костей, а также в изменениях фи- 
зико-механических их характеристик, прежде всего упругости, вязкости 
и пластичности. Такие изменения лежат в основе процессов онто- и 
филогенетического формирования конечностей четвероногих различных 
морфо-функциональных типов, и поэтому исследование их представляет 
‘большой интерес для понимания путей и закономерностей становления 
этих типов. В последние годы необходимость таких исследований увели- 
чивается в связи с развитием бионики. Речь идет о физическом модели- 
ровании принципов строения и функционирования конечностей в целом 
и их отдельных механизмов. Несомненно, принцип «выращивания» дета- 
лей механизмов и машин, обладающих при минимальной массе опти- 
мальной выносливостью к определенному комплексу нагрузок, будет 
позаимствован у живой природы и позволит значительно (в два-три 
раза) облегчать механизмы и машины без ущерба для их прочности. Но 
этому должен предшествовать не только качественный, а и количествен- 
ный анализ структуры и физико-механических характеристик биологи- 
ческих прототипов. 

В работах ряда исследователей (Гильзен, 1896; Ѕтіїіһ, \а|тз1еу, 
1957, 1958; $е4Ип, 1965 и др.) установлено, что на кости конечностей при 
движении животного действуют почти мгновенные силы и что вещество, 
из которого построены кости, под действием этих сил ведет себя как 
упруго-вязкий материал, т. е. сопротивление костного вещества дей- 
ствующим на него силам зависит от скорости деформирования. Реология 
(наука о деформировании и текучести материалов) доказывает, что нель- 
зя получить достоверное представление о способности упруго-вязкого ма- 
териала сопротивляться быстро меняющимся во времени нагрузкам, ес- 
ли не учитывать наряду с упругими, также их вязкие свойства (Рейнер, 
1965). Поэтому при исследовании реакции кости на динамические на- 
грузки нужно глубоко изучать ее упруго-вязкие свойства. Однако таких 
исследований проведено еще мало. Одна из причин — слабая разработка 
их методики. Существующие же методики определения упруго-вязких 
свойств технических материалов по разным причинам не могут быть не- 
посредственно использованы при биологических и бионических исследо- 
ваниях. 

В исследованиях лаборатории бионики отдела эволюционной морфо- 
логии Института зоологии АН УССР изучению упруго-вязких свойств 
костей отводится большое место, и поэтому мы задались целью разра- 
ботать простую методику исследования данных свойств костной ткани 
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и показать пути анализа реакции костного вещества на воздействие сил, 
характер изменения которых во времени сложный, а сами эти силы в пер- 
вом приближении могут быть сравнимы с силами, действующими при 
функционировании конечностей. 


Постановка задачи 


Некоторые исследователи (Ѕейііп, 1965; Аѕсеп2і апа Вописсі, 1967, 
1968) показали, что реакцию =(#) костного вещества в пределах упругих 
деформаций на силы 0(#), действующие на него, можно изучать с по- 


Ёк 


у &0 
ГА 


Рис. 1. Линейная упруго-вязкая мо- Рис. 2. Многоэлементная линейная упруго-вяз- 
дель костной ткани (объяснение в кая модель: 

тексте). Е, Е, Е, Ез...Е „— коэффициенты упругости; Ти ТП... 

Е Рн коэффициенты вязкости; 0 — действующая сила. 


мощью линейной упруго-вязкой модели, состоящей из тела Кельвина с 
параметрами Е, и п, и соединенного с ним последовательно элементом 
Гука с упругой константой Е, (рис. 1). Ставится задача: определить 
упругие и вязкие параметры этой модели и произвести анализ реакции 
ее на воздействие динамических сил. 


Решение 


Комплексный модуль 


Пусть к линейной модели, описывающей упруго-вязкое поведение 
'материала (рис. 2), приложена синусоидально изменяющаяся сила 
6 = Ш ®ё с угловой частотой Ф. По истечении времени, достаточ- 
ного для того, чтобы можно было пренебречь начальными условиями, 
деформация, или реакция модели, будет меняться по тому же закону 


Рис. 3. Графическое изображение ре- 
акции линейной упруго-вязкой моде- 
ли на действие периодической силы: 


1 — реакция модели == 5ѕіП(оѓ — Ө); 2— 

сила с периодом Т и частотой 0; Ө — фа- 

зовый сдвиг между силой и реакцией мо- 

дели на нее; ©; и ё — соответственно наи- 

большие амплитуды силы и реакции моде- 
ли при ®=соп$(. 


Амплитуды силы и деформации 


с угловой частотой о и запаздывать по фазе на угол Ө (рис. 3) (Кар- 
ман и Био, 1948; Альфрей, 1952). В таком случае сила и деформа- 
ция окажутся связанными общим комплексным числом, или т. н. 
комплексным модулем У (іх) * (Бленд, 1965). Для модели (рис. 2) 


* Символ ({ «) указывает на то, что данная функция комплексная и зависит от 
параметра о. 
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это число запишется как 


1 
ү 22 8 (1%). 38 : 1 
бо = = (іо) а ( А 
1 у 1 
Е и 
Е, Ер + ть 
&=1 

Абсолютная величина его выражается соотношением 

| (25) |= ИМ = №, (2) 


а фазовый сдвиг между силой и деформацией — углом 
М 
Ө = агс4е — , 3 
Еф (3) 


где М — действительная и № — мнимая части комплексного модуля (1); 
они являются функциями упругих Е, Ё, Ёз,... Е, И вязких 11, 1, ... 12 
параметров модели на рис. 2, а также частоты о. 

Из уравнения (1) и рис. 2 видно, что элементы Ёо и Е, комплексно- 
го модуля Ү(і~) определяют упругость материала, в то время как эле- 
менты п, определяют его вязкость. Следовательно, при ® =сопѕі, более 
упругий материал имеет большие величины Ро и Ё, и меньшие величи- 
ны п, и наоборот. Угол Ө, определяемый формулой (3), отражает вели- 
чину энергии рассеяния. В технике он называется механическим углом 
потерь, а само явление диссипативного рассеяния энергии материалом 
принято называть внутренним амортизатором (Писаренко, 1955; Тимо- 
шенко, 1967). Записав формулу (1) в виде 


=, у 9 (4) 

а) = (4) 
легко увидеть, что простое деление амплитуды заданной силы о (іо) 
на абсолютную величину модуля У ({®) дает амплитуду реакции мате- 
;Риала с (і є). 

Из изложенного выше вытекает важный факт: если для данной мо- 
дели комплексный модуль Ү (ію) известен, то реакция материала на 
синусоидальные напряжения данной частоты вполне определима. Ее ам- 
плитуда и фазовый угол вычисляются по формулам (3) и (4). 


Преобразование Фурье 


Ниже будет показано, что для получения аналитического выраже- 
ния комплексного модуля Ү (іо) понадобится формула преобразования 
функции времени =({) и с (і) в частотные функции о (іо) и о (іє). 

Известно, что любая функция времени | ({), подчиняющаяся опреде- 
ленным условиям, может быть преобразована в функцию $ (і о) от не- 
которого параметра ®, причем между обеими функциями существует 
следующая связь: 


< 


еба) = 2 \ Ре аё. (5) 


Частотная функция Ф (і о) называется преобразованием Фурье функции 
времени {(Ё). Так как в дальнейшем изложении мы будем оперировать 
функцией времени, не принимающей отрицательных значений, то нами 
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используется выражение (5) в виде 
е 0) = а | 0 еа (6) 
2п 
0 
Пользуясь соотношением (6), преобразуем производную функцию вре- 
мени $’({) в частотную $'(2) 


бы) = | Г-н (7) 
2= . 
б 
Интегрируя по частям выражение (7), имеем 


е а) = 100) етер + | 10) еа. (8) 
7) 
Чтобы соотношение (8) имело смысл, необходимо соблюсти условие 
[0 е] 0 при #- ос. 
Тогда, учитывая соотношение (6), получаем 


(бы) = —— в (а) + 0), (9) 


где (0) — значение функции #(/) при #=0. 

Уравнения (6) и (9) являются теми искомыми, которые необходимы 
для отыскания функций 2 (і о) и =(2 о). В несколько измененном виде 
‚они известны в литературе (Альфрей, 1952). 


Определение комплексного модуля, 
упругих и вязких параметров механической 
модели 


А Поскольку комплексный модуль позволяет определять реакцию 
упруго-вязкого материала на переменные силы, а также анализировать 
‚упругие и вязкие характеристики материала, необходимо уметь его вы- 
ЧИСЛЯТЬ. 

Существует два основных метода получения комплексного модуля 
материала для данного типа нагрузки: один основан на серии динами- 
ческих тестов, а другой— на использовании простых статических опы- 
тов (Ѕітһа, 1943; Альфрей, 1952). С помощью динамических тестов изу- 
чают реакцию материала на ряд простых синусоидальных напряжений 
при варьировании частоты о от очень низких до очень высоких величин 
и комплексный модуль получают как дискретную пару величин числа 
и фазы. Этот метод позволяет графически изобразить два компонента 
комплексного модуля, но не дает аналитического выражения его. Поль- 
зуясь вторым методом, заранее известное °({) (его преобразование 


с(і о) по формуле (6) имеется) прикладывается к испытываемому мате- 
риалу и регистрируется его реакция =({). Если есть возможность преоб- 
7с (іо) 
= (ѓо) 
‚дает комплексный модуль в аналитической форме в функции частоты о. 
Таким образом, очень простой опыт может дать полную информацию о 
комплексном модуле. Единственная проблема, возникающая при вычис- 
‹лении комплексного модуля таким путем, заключается в получении пре- 


разовать є (і о) полученного ответа &({), тогда отношение прямо 
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образования Фурье по формуле (6) зарегистрированного ответа =(#). 
Она может быть решена с помощью вычислительных машин. Если же 
по каким-либо причинам невозможно пользоваться вычислительной ма- 
шиной, то используют известные графо-аналитические методы апрокси- 
мации кривых (Анго, 1967). Большинство из них основывается на 
апроксимации ответа =(!) с помощью ряда простых функций, преобра- 


5 
Еч 
35 ——— 
Ё А 
За ТО 
& 
ЕВ 
9 
ее ГРА 
С. 
Кк 
ь Время, сек 
а $ Не 
СА 
Рис. 4. Статическое испытание костной ткани Рис. 5. Электрический ана- 
на растяжение: лог механической модели, 
а — схема опыта; б — графическое изображение ре- изображенной на рис. | (объ- 
зультата; / — испытываемый образец; 2 — микроэк- яснение в тексте). 


стензометр; 3 — усилитель сигнала; 4 — регистрирую- 

щее устройство; о, — растягивающая сила; е — мгно- 

венная реакция костной ткани на силу 01; ё — полная 
реакция (деформация) костной ткани. 


зования = (і ®) которых существуют, по формуле (6). В таком случае 
общее преобразование отыскивается как сумма преобразований от про- 
‚стых функций. 

Предположим, испытание материала проведено согласно схеме 
(рис. 4а). и получена реакция его на напряжение а в виде графика 
(рис. 4 6). После графического дифференцирования функции :(#) апрок- 
симируем кривую =’(!) рядом экспоненциальных кривых. Операция 
апроксимирования выполняется таким образом: 1) вычерчиваем кривую 
#корости деформации в полулогарифмической системе координат {— Ш =’ 
(скорость деформации откладываем на вертикальной логарифмической 
шкале) и нижний участок кривой апроксимируем прямой линией; 2) раз- 
ность в ординатах между кривой и прямой линиями подсчитываем как 
функцию времени и результат вычерчиваем в упомянутой выше системе 
координат; 3) нижний участок вновь полученной кривой скорости дефор- 
мации апроксимируем прямой линией и вторично подсчитываем разность 
в ординатах между кривой и прямой линиями как функцию времени и 
результат вычерчиваем в виде графика (см. п. 1); 4) переходы операции 
апроксимирования повторяем до тех пор, пока остаток разности между 
кривой и прямой линиями можно полностью представить прямой лини- 
ей; 5) каждую из прямых линий (экспонент) определяем формулой типа 
А,е— * . Сложив уравнения всех прямых, получаем аналитическое вы- 
ражение для скорости деформации в форме 


60) = Ў Де. (10) 
1—1 


Значения А; находим по пересечению прямых с осью ординат, а вели- 
чины Ё; — по наклонам прямых. 

Для получения преобразования Фурье функции = ”({) используем 
формулу (6) 
5—1869 
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Е” ыы [| (Але-& + Азе-ы - ... + Ае ) еі Е = 
6 
== | [А енн А, ене 4... А, АУ | и 
п 
0 
ТА А, Ал. 1 
[а т 05 ` Е в 
а по формуле (9) получаем преобразование с (і ~) деформации =(й). 
Ен ие. А, 
9 а Ея № іо а» ё + і 2) +5]. (12) 
Если о (4) = о, для Ё>0и 
ы с (0) = 0 для # < 0, тогда з (№) = 1 . (13) 
9] 
Используя основное уравнение (1), получаем 
аа рен 
Ү (іо) ас д. х УЕ . (14) 
п \ + Ее Е с +“ 


Определение упруго-вязких констант материала легко выполнить, если 
уравнение (14) записать в виде 


1 

фо) — и, (10 
Ү(іч) = ? г" ) 

22; 2паџ: А Эк 1. пар `` 2 о д а. 

р іь 0 5 і ———— Ию — 

„А, А; А, А, Ал Ал о 

а уравнение для комплексного модуля модели на рис. 2 — в форме 
: 1 
Ү (іо) = — (16 
(о) Е Е еа 
Е, + №1: Ез + ют = Еһ + іәт Ё, 


и сравнить их. Нетрудно видеть, что соответствующие действительные и 
мнимые части равны. Так, 


Е, = 23а Е, = 1. а Вара. Е ‚ = 
1 з п 0 
: (17) 
е 2с РО 221 РЕ 2па; 
РА а, В Я 


Итак, с помощью методов математики и простых опытов, сущность 
которых изложена выше, можно вычислить комплексный модуль (1), 
а также вязкие и упругие параметры модели, описывающей упруго-вяз- 
кие свойства материала (17). 


Пути анализа реакции материала 
на динамические силы 


Анализ реакции модели, описывающей упруго-вязкие свойства дан- 
ного материала, на динамические силы можно выполнить, воспользовав- 
шись аналитическим выражением (4). Но такой путь применим, когда 
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элементов, из которых состоит модель, мало (три И меньше). Анализ 
реакции многоэлементных моделей удобнее выполнить на электронной. 
аналоговой машине или, если по каким-либо причинам такой машиной 
воспользоваться нельзя, на электронном приборе, построенном специаль- 
но для этих целей. Располагая данными о структуре, упругих и вязких 
параметрах механической модели, такой прибор нетрудно построить, 

1(0) костной ткани в пределах упругих деформаций на 
динамические силы: 


1 2 3 1 — входная функция; 2 — оператор; 3 — выходная функ- 
ция. 


Рис. 6. Блок-схема прибора для исследования 


если уметь составлять схему электрического аналога механической мо- 
дели. Основные принципы составления таких схем изложены ниже. 

Известно, что каждому простому типу механической деформации 
соответствует простая электрическая цепь, обладающая теми же свойст- 
вами (Карман и Био, 1948). 


Электрические аналоги 


Механическая система Электрическая система 


Напряжение (1) Напряжение У(2) 
Деформация (і) _ Заряд 4(1) 
Скорость деформации _ Сила тока КЁ) 
Е ’(#) 
1 
Податливость 1= Е Емкость С 
Вязкость 1 Сопротивление К 


Напряжения складывают- 
ся при параллельном сое- 
динении 

Скорости деформации 
складываются при после- 


Напряжения складывают- 
ся при последовательном 
соединении 

Токи складываются при 
параллельном соединении 


довательном соединении 


Используя эти соответствия, механическую модель, описывающую 
упруго-вязкие свойства костной ткани (рис. 1), можно представить. 
в виде электрического аналога (рис. 5). Такая модель может быть реа- 
лизована в приборе, блок-схема которого представлена на рис. 6. При- 
бор позволит произвести исследование реакции костной ткани на дина- 
мические силы в условиях, которые в лаборатории создать трудно или 
почти невозможно. 


Заключение 


Описанный нами метод исследования упруго-вязких свойств трубча- 
тых костей прост, доступен, не требует дорогостоящих приборов и дает: 
достоверные результаты. Кроме того, он универсален: при определенных. 
условиях позволяет исследовать не только костную ткань, но и хрящи,. 
связки и другие биологические ткани. 
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